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1. Introduction 



There has been much speculation that stable topological defects such as monopoles, 
strings, or domain walls may have played an important role in the early universe. Unfortu- 
nately, such defects only arise in extensions of the standard model which are rather poorly 
constrained by experiment. It was discovered recently that soliton-like string solutions 
can exist in an unphysical version of the Weinberg- Salam model with the SU(2)l gauge 
coupling set to zero fl|] , || . Although strictly speaking these solutions have no topological 
charge (the vacuum manifold is simply connected), they may still be characterized by a 
"winding number" N which characterizes the variation of the phase of the Higgs field as 
one encircles the string. This winding is associated with the fact that the gauge orbit space 
is not simply connected. 

Semi-local strings also provide solutions to the physical version of the model with 
non-zero gauge coupling [[|. String solutions in the Weinberg- Salam model had also been 
discussed earlier in |4]]. For physical values of the coupling constants the solutions with 
N = 1 are unstable ||. The linear stability analysis of || does not indicate whether the 
instability leads to a dissolution of the string, or whether it leads to a modified solution of 
lower energy per unit length. In || it was shown that the instability can be viewed as due 
to a condensate of charged vector bosons in the core of the string, a phenomenon which 
has been investigated in the context of electroweak symmetry restoration in magnetic fields 
. In some qualitative arguments were given suggesting that a modified solution may 
exist when sin 2 6w = 0. As we will show, the perturbation which seems most likely to 
destroy the string has improved stability for winding N > 1, suggesting that higher winding 
strings are a natural place to look for a modified solution involving W condensation. In this 
paper we will show that for A > 1 there is an instability due to W condensation which leads 
to an apparently new string solution which we term the "Cinderella String" at the physical 
value sin 2 6w ~ 0.23 when the Higgs mass is less than the Z mass and we will construct 
this solution numerically. We will then explicitly construct a gauge transformation that 
converts this Cinderella string into an ordinary electroweak string with winding number 
one less. 
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2. Semi-local and Electroweak strings 

The action for the bosonic fields in the Weinberg-Salam model is 

S = J d 4 x (L> M $)t(D M $) -Ip^p^ - ^G a ^G a ^ -A($t$_ ^2/2)2 



(2.1) 



where P M $ = (V fl -igW fJ a a a /2-ig' 'B^/2)$. and W a (a = 1, 2, 3) are gauge potentials 
for U(l) and SU(2)l respectively, with field strengths = d^B u — d u B^ and G a v = 
9/iW u a — d v W^ a + ge ahc W b W u c . Our conventions are such that e 123 = 1; a a are the 
Pauli matrices. Also, W± = {W^ 1 ± i W 2 )/V2 , = cos9 w W* - sin^S^ , A M = 
sin^W 3 + cos 9 w B fl where sin 2 9 W = g' 2 /a 2 rs 0.23 with a 2 = g 2 + g' 2 . In what follows 
we will work in rescaled cylindrical coordinates (t,z,r,9), where r = y/Xrjp (p being the 
usual radial coordinate), and use a coordinate basis rather than the orthonormal basis 
often used. 

The equations of motion which follow from ( |2.1|) are 

2 

D^D^ + 2A($ f $ - ^-)$ = 

V fM F t " / + —(&D v ®-(D v ®y®) = (2-2) 
V^G^ a + gt abc W^G^ c + —{&<r a D v $ - (^f)V a $) = 

It was shown in for 9 W = n/2 and in [[| for arbitrary 9 W that these equations admit a so- 
lution describing a "Z-vortex": Z = 2N(P(r) - l)/a, $ = {<f> w fa) = rj(0, X(r)e lN0 )/V2 
(all other fields set to zero), where X(r) and P(r) are the solutions to the Nielsen-Olesen 
equations (' = d/dr): 



-X"-- + 1 -^ Y ^+X(X 2 -1) = 

P> 2 (2 - 3) 
- P" + — + -PX 2 = 

r p 

satisfying the boundary conditions X(0) = 0, P(0) = 1, and -X"(r) — > 1, P(r) — > 0, as 
r — > oo. The parameter /3 is the ratio of the square of the scalar and vector masses, 

P = Wffiggs/™! = 8A/a 2 . 

Note that this is not a "topological" string in the usual sense. The vacuum manifold 
of the Weinberg-Salam model is S 3 and it has no non-contractible loops. However it has 



been shown |l],0,|3|,[5|] that, for certain ranges of values of the Higgs mass and the Weinberg 
angle, these configurations are stable to small perturbations. 

A special case where stability is easy to understand is N = 1, cos9 w = 0. In that 
case the SU(2) gauge fields decouple (g = 0) and we have what are called semi-local 
strings, with symmetry SU(2) g i b a i x U(l)i oca i and the same vacuum manifold as before. 
Unlike in the case of Nielsen-Olesen strings, where iV = 1 is automatically stable to small 
perturbations, stability here is decided by a competition between gradient and potential 
energy: when (3 > 1 there are nearby configurations with a non-zero <p u at the core whose 
increased gradient energy is more than compensated by a lower potential energy (from 
the last term in the action) [pj. For (3 < 1, however, semi-local strings are perturbatively 
stable |I],^,I0]. Physically this is because the variation of the Higgs field at infinity in 
the S 3 vacuum manifold lies on the orbit of the U(l)i oca i gauge symmetry. To unwind 
the string must move off this orbit, but this raises the energy because the energy due to 
the angular dependence of the Higgs field can no longer be canceled by a U(l)i oca i gauge 
field. When cos (9^ 7^ the barrier to decay becomes finite or may disappear and one must 
investigate the stability numerically. 



3. W condensation and large-iV vortices 

There are two distinct physical mechanisms that can lower the energy of an electroweak 
string. The first one we have just discussed: an instability occurs if it is energetically 
preferable for <j) u to develop a non-zero value at the origin. The second one, pointed out 
in is the formation of a condensate at the core of the string. For N = 1 there is 
in fact no clean separation of the known instability into these two types of perturbation. 
The equations of motion require a non-zero <p u at the origin as soon as there is a non-zero 
W condensate. For N > 1 however there is a clear distinction between the two. The 
first instability seems likely to destroy the string, as in the semi-local case, since once the 
Higgs field moves away from zero at the origin there is no reason for it not to move back 
to its vacuum value everywhere. On the other hand one might expect the second form of 
instability to lead to a new string solution since it does not seem to change the winding 
structure in the Higgs field. It was even suggested in |§ that a "VF-condensed" solution 
might be more stable than the electroweak string. In what follows, we investigate these 
claims. 
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It can be shown that, inside the core of a large- N vortex, the Nielsen- Olesen equations 
( pT3l) reduce at leading order in 1/N to 

P' 

-P" + — = 

\, p (3-1) 

£ r 

with solution 

P = 1 - pr 2 ; £ = Cre~ pr2/2 (3.2) 

where £ = X 1 /^. (The constant C can be estimated by demanding that £ ~ 1 when P pa 
which gives C ~ ^/pe; its precise value will not be important in what follows) . In deriving 
these expressions we have used the fact that < £ < 1 and kept leading terms in N. We 
now verify that the £' and £" terms dropped are indeed subleading by calculating the N 
dependence of the coefficient p explicitly. 

Since X pa (per 2 ) N ^ 2 e~ Npr I 2 in the core, transition from X = to X = 1 occurs 
fairly abruptly at r pa 1/ ^fp with a width of ~ 1/ y/Np. To study the transition region we 
introduce the variables 



u=y/Wp{r-^-) ] 7 = ^; Q = ^- (3-3) 



The leading order behavior of the equations (|2.3|) then takes the form 

-X" + AXQ 2 + ^ 2 X(X 2 - 1) = 



Q + -irQ x = 



with 

O = I 

w>0 1 1 « > 

The equations (EO]) have a first integral 



Q = {" U ' W<< ° X = {°' W<< ° (3.5) 



-i(X') 2 - 4> (Q') 2 + 2^ 2 Q 2 + y(^ 4 " 2X 2 ) = Const (3.6) 

The boundary conditions at u = +oo fix the constant to be — 7 2 /4, while the boundary 
conditions at u = — oo give 



We thus find p = l/{2Ny/]3) when A" >> 1. This shows that we were justified in dropping 
£' and £" terms to derive ( |3.2| ). It also implies the large- N scaling of the core radius 
r core ~ \/2N(3^ (the width of the transition region stays of order unity). 

Now consider a (z- and 9- independent) perturbation in <p u which is non-vanishing at 
the origin 



* _ V_ ( /(r)e°* \ 
y/2\X{r)e im ) 



/(r)~/ + O(r 2 ) asr^O (3.8) 



The corresponding W perturbations decouple if f(r) falls to zero while inside the core of 
the string, (since they always appear multiplied by X « in the / equation). Inserting 
the explicit form for P yields 



t „ f (2cos 2 ^-l) 2 2 fQ 2 
r 4/3 r 1 \Xr] 



+ V J; ^ r*f+(^-l)f = 0. (3.9) 



The equation can be rescaled to depend on only one parameter by the change of variables 
r 2 = 2y/]3r 2 /\2 cos 2 6 W — 1|, giving 

-/"(r)- ^P + (f 2 -A)/(f) = 0, where A = 2 ^f - ( ^ - l) (3.10) 
r \2 cos z V w — 1\ \Ar] z J 

The solution is f(f) = foe~ r I 2 and it corresponds to A = —2 (which is the lowest 
eigenvalue since f has no nodes). In principle there could be lower eigenvalues coming from 
the other equations (those of the W perturbations) but one can show that this is not the 
case. The line of critical stability is Q = 0, or 

]3 = \2 cos 2 6 W - 1\ (3.11) 



(See Figure 1). This argument shows that the stability of bare electroweak strings with 
respect to this particular perturbation is better for large- strings than for their N = 1 
counterparts (for sin 2 9 W = 0.23 we find the critical line at mniggs ~ 4:8.6GeV). Note that 
consistency requires that / reach zero before X becomes non-zero; for the physical 9 W this 
is true to within 1% when A^ > 23. 
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Fig. 1: Stability of electroweak strings to perturbations with a non-zero <p u at 
the center. The dashed line shows critical stability for large- iV strings to the 
specific perturbation of eq. (3.8). The solid line shows critical stability (against all 
perturbations) for the N = 1 string, in agreement with H. 



We now consider stability with respect to W condensation. Recall that inside the 
string core P only deviates from quadratic behavior when X is substantially different from 
zero. Thus, for large N there is a Z magnetic field which is approximately constant 
throughout the core with a magnitude given by ANp/a = 2/ay/]3. It is known from the 
work in in a different context that such a uniform magnetic field is unstable to the 
formation of a W condensate. The analysis in ]7| shows that such an instability exists if 
the magnetic field is greater than a critical value determined by the value of the W mass. 
Since the W fields are essentially massless in the core of the vortex at large N, the field in 
the core is large enough to cause this instability. The static energy is minimized by having 
a condensate with a constant W magnetic moment in the core. In cylindrical coordinates 
this suggests an ansatz of the form 

W± = ±iV\rjv(r)e iv0 ; W± = m{r)e lv \ (3.12) 

where v is a positive integer. An extension of the analysis in || for arbitrary N shows that 
an efficient way to lower the energy which is compatible with the equations of motion is to 
have v = 1 and the leading r behavior given by v(r) = vq + 0(r 2 ), m(r) = —v$r + 0(r 3 ). 


1 Note that the angular dependence of the condensate differs by a sign from || and the radial 
dependence is one power of r less than that given in || for N = 1. 
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Moreover, analysis of the equations of motion with this angular dependence shows 
that it is not consistent to set the upper component of the Higgs field, 4> u , to zero. An 
upper component of the form y(r)e l( ' Ar_1 * )6 ' is required, with Y(r) = yor 1 ^ -1 + 0(r N+1 ). i 

We thus see that at large iV the string has improved stability with respect to pertur- 
bations of 4> u and has an instability due to W condensation. In the following section we 
follow numerically the instability due to W condensation to find a modified string solution. 



4. A W dressed solution 

Following the previous discussion we consider the following ansatz for a W condensed 
solution 

d> ee fa, ft) = \ (Y(r)e^ N -^,X(r)e iNd ) 

(4.1) 



W r = iV\r]v(r)e id ; W e = m(r)e l9 

Z e = — (P(r) - 1); A e = —A(r) 
a a 



With this ansatz the energy is given by 



E=\]r d r\(Y>- ^ 



9vX^ + {x , + gvY )2+ i_ {x2 + Y2 _ i)2 



V2 



+ \(NPX- 



gmY 2 (3a 



P r 2 



+ 



w! + v ( 1 + %N(P - 1) + ^-NA 
or cr 



- ~^[g 2 NP' + gg'NA']vm + 







NP'Y fNA'V | g 2 a 2 v 2 m 2 



r 



r 



(N -l)Y-( g2 -f\ N(P - 1)Y - %NAY - ^£ 



-i 2 



(4.2) 



To find a solution of the equations of motion (|2.2| ) we used a variational method to 

look for minima for the energy (|4.2|) given the ansatz ( |4.1| ). Each of the fields X l / N , 
yi/(Ar-i) 5 p ; ^ m and 

v in the energy was written as a finite Fourier series running from 
r = to a sufficiently large outer radius r out (typically twice the core radius), and the 
harmonics' coefficients varied until the energy was minimized. The harmonic series were 
constrained to satisfy the boundary conditions X x / N ~ r, yVC-W-i) r^r,Ar^r 2 ,m^r, 
v ~ const, P — 1 ~ r 2 at r = 0, and P was forced to be zero at radius r out . The minimum 

2 For general is < N we find Y ~ e KN-u)e r N- v as r ^ . 
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energy configuration was then refined by increasing the number of harmonics, until the 
energy stayed constant to one part in 10 5 . We then substituted these configurations into 



the equations of motion (|2.2|) to confirm that we had found a solution. Working with 15 
basis functions for each function and using 60 discrete radial points in the integrations, we 
found that the equations of motion were satisfied except for high-frequency residuals. In 
figure 2 we compare, as a function of (3, the energies per unit length of the VF-condensed 
solutions with N = 2, 4 and of bare electroweak strings with N = 1,2, 3, 4. 



Fig. 2: The energy per unit length for the electroweak strings (solid lines) as a 
function of (3 and N, and for the Cinderella strings with N = 2 and 4 (dots). 



We have also looked for solutions that minimize the energy for TV = 1. For the range 
of (3 and sin 2 9 W shown in figure 1 we found that the bare electroweak string was a solution 
and was stable to the development of a W condensate. Our results are in agreement with 
those of . For values of (3 and sin 2 9 W outside this range the program would only converge 
to a vacuum solution, consistent with the lack of existence of a non-trivial W condensed 
electroweak string with N = 1. 



5. The clock strikes twelve 

Inspection of the energy per unit length of the Cinderella string with winding N and 
the bare electroweak string with winding N — 1 shows that they are the same to our 
numerical accuracy. This strongly suggests that they are in fact the same solution up to a 
gauge transformation. To show that this is the case we construct a gauge transformation 
relating the two solutions and verify numerically that they are in fact equivalent. To 
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construct the gauge transformation it is simplest to first focus on the Higgs field. For the 
Cinderella string it takes the form given by fl4.1| ) with asymptotic behavior 



X(r)->1, y(r)->0, r^oo 

(5.1) 

X(r) -> xor^, y(r) -> j/o^ -1 , r -> 0. 

This may be put in the form of the Higgs field for a bare electroweak string of winding 
N — 1 , $ = G $ by the gauge transformation 

G= fcosipe ie -sinV> \ / g 2 n 

y sin-0 cosipe~ %e J 

where 

cosibjr) = , K J = (5.3) 
^JX 2 (r) + Y 2 (r) 

We have verified numerically that the resulting Higgs field profile is the same as that for a 
bare electroweak string of winding N—l. We have also checked that the gauge transformed 
gauge fields, 

W r = iV2 (v(r) - ^Vj e~ l ° 

W e = V2( m(r) cos 20 + (W$(r) + e -« (5-4) 

V 9 V2 / 

2 

W/f = W|(r) cos 2-0 - V2m(r) sin 2*0 + - cos 2 -0, 

are those of an electroweak string with winding N — l, that is, VF r = VFe = A# = and 
Zq = 2(N — l)(P(r) — l)/a with P(r) the Nielsen-Olesen solution for winding N — l. 



6. Conclusions 

We have shown that for the physical value of sin 2 9 W electroweak strings with winding 
N > 1 have an instability associated with the formation of a condensate of vector 
bosons in the core of the string. By minimizing the energy with an ansatz compatible 
with this condensate we were able to find a W condensed Cinderella string. However this 
new Cinderella string is an illusion, it is gauge equivalent to a bare electroweak string with 
winding N—l. 

By following the W condensate instability numerically to a new solution we have 
shown that there are no new W condensed electroweak string solutions within the ansatz 
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we have used. While we cannot completely rule out the possibility of more complicated 
stable string solutions in the Weinberg- Salam model, our results make it quite likely that 
none exist. This leaves only the N = 1 electroweak string which, as has been shown in a 
previous analysis || and reconfirmed here, is not stable for physical values of the Weinberg 
angle. 
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